Double poles in Lattice QCD with mixed actions 
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ABSTRACT 

We consider effects resulting from the use of different dis- 
cretizations for tlie valence and the sea quarks, considering Wil- 
son and/or Ginsparg- Wilson fermions. We assume that such ef- 
fects appear through scaling violations that can be studied using 
effective-lagrangian techniques. We show that a double pole is 
present in flavor-neutral Goldstone meson propagators, even if 
the flavor non-diagonal Goldstone mesons made out of valence or 
sea quark have equal masses. We then consider some observables 
known to be anomalously sensitive to the presence of a double 
pole. We find that the double-pole enhanced scaling violations 
may turn out to be rather small in practice. 
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1 Introduction 



The use of mixed actions in Lattice QCD (LQCD) is at present gaining 
popularity. With "mixed action" we refer to any simulation in which different 
lattice fermions are used in the valence and sea sectors. (If only the bare 
masses differ, we are dealing with a partially-quenched (PQ) theory pp.) 
The main reason is that the generation of dynamical-fermion configurations 
is very expensive, while for many quantities the "good quality" of symmetries 
(such as chiral and flavor or taste symmetries) is more important in the 
valence sector than in the sea sector. An example of this are the power- 
like subtractions needed for K ^ it, for which good chiral symmetry in the 
valence sector is essential, but not that in the sea sector. 

However, the use of mixed actions raises field-theoretical questions. At 
non- vanishing lattice spacing a a mixed-action theory is not unitary, and the 
question is what happens to unitarity in the continuum limit, as well as how 
big the effects of the violation of unitarity are at the lattice spacings used in 
actual simulations. We note that similar issues arise for many discretizations 
of LQCD which do not employ mixed actions, such as improved actions and 
actions with Ginsparg-Wilson (GW) fermions. 

Here we address the second question, the size of effects at non-zero lattice 
spacing. We extend the notion of universality to assume that unphysical 
effects due to the use of a mixed action vanish in the continuum limit, and 
are controlled by positive powers of a. We further assume that effective-field 
theory (EFT) techniques can be used to investigate the behavior of mixed- 
action theories at low energy. 

The most salient phenomenon which occurs if the valence quarks don't 
match the sea quarks is the appearance of a double pole in flavor-neutral 
Goldstone meson propagators, with a residue R proportional to a^, as will 
be explained below. If also the renormalized quark masses are not equal, 
there is also a contribution from the quark mass difference, R oc cia^ -|- 
C2{'>TT'sea — ^valence) |I] • Slucc the double polc represents a dramatic change 
of the infrared (IR) behavior of the theory, it is important to investigate the 
most serious consequences of the double pole. As an example of this, we 
consider the enhanced finite-volume effects in J = pion scattering. For a 
more detailed account of this work, see Ref. Pj. In our work, we deal only 
with the cases in which the fermions in both sectors are Wilson-like or GW 
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fermions. 



2 The double pole 

According to our assumptions, a mixed-action theory is a PQ theory in the 
continuum hmit |2j. We therefore start with continuum chiral perturbation 
theory (ChPT) at lowest order, described by the PQ lagrangian |1| 

Leant = If str(a^Ea^Et) - IfBo str(SMt + MSt) , (1) 
o 4 

where S = exp(2i$//) is the non-hnear meson field, / and Bq are 
the well-known lowest-order low-energy constants (LECs), and M = 
diag(mt,, m^, . . . , m^, m^, . . . , m„, rriy, . . .) the quark mass matrix, with K va- 
lence quarks of mass rriv, N sea quarks of mass nis, and K ghost quarks of 
mass rriv'^ The relevant chiral symmetry group is thus SU{K + N\K)l x 
SUiK + N\K)n. 

The next step we need is provided by the systematic analysis of Ref. j3] , 
where mixed actions with a Wilson sea and GW valence quarks were con- 
sidered in ChPT up to order a^. It is very easy to generalize that work to 
the slightly more general case where the valence and sea sectors can be any 
combination of Wilson-like and GW fermions.^ 

The first observation is that at a 7^ the symmetry group is smaller. For 
GW valence quarks and a Wilson-like sea the symmetry group is 

[SU{K\K)l X SU{K\K)r] X SU{N) (GW - Wilson) , (2) 

while the case of Wilson-like valence and sea quarks (but with unequal lattice 
actions) the symmetry group is even smaller: 

SU{K\K) X SU{N) (Wilson - Wilson) , (3) 
^For mixed actions with staggered sea, see Ref. 

^We take degenerate valence and degenerate sea quarks, but it is straightforward to 
generalize. 

•^With "Wilson-like" we refer to any type of Wilson fermion, including for instance 
clover and twisted-mass Wilson fermions. Our conclusions can be straightforwardly gen- 
eralized to these situations. With GW fermions we refer to overlap fermions or domain- wall 
fermions with a small enough residual mass. A potentially interesting special case is the 
combination of untwisted Wilson sea quarks and twisted Wilson valence quarks. If this is 
done by only twisting the mass matrix of the valence quarks, this is a PQ theory, and not 
a (genuinely) mixed theory. The relevant symmetry group in that case is SU{K -f N\K). 
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because all chiral symmetries are now broken. 

Up to order a} there are a number of new operators which appear in the 
effective lagrangian, and we only give the ones here which are relevant for 
our investigation, 

^Lw = -^(7..(str(P,(S-St)))2 + 7,,(str(P,(S-St)))2 (4) 

+7,,str(P,(S - St))str(P,(S - St))) , 

in which '^vv,vs,ss are new LECs, and P„_s are projectors on the valence and sea 
sectors, respectively. We note that ■jm = Jvs = when the valence quarks are 
GW, because these operators break chiral symmetry, and thus are excluded 
by the symmetry group (jS]), but not by the group 0. 

In order to obtain the meson propagators, we need to expand the la- 
grangian to quadratic order in the field $. If we integrate out the sea-?]' 
(which is heavy because of the singlet part of the t]' mass), we have that 
str(<l>) = str((P„ + Ps)$) = 0, and Eq. (jU reduces to quadratic order to 

6Lw = (7™ + 7ss - 27„,) (str(P„$))2 + . . . . (5) 

Recognizing that str(P^$) is nothing else than the "super-?]'" field of 
quenched QCD, it is obvious that 6Lw leads to a double pole in flavor- neutral 
propagators. 

In the flavor non-diagonal sector, we find (tree-level) Goldstone meson 
masses jl] 

Ml = 2Bo,m, + 2Wo,a + 2(3,a^ , (6) 
Ml = 2Bosms + 2Wosa + 2/3,0^ , 

where we now allowed for the fact that with the reduced symmetry groups 
the LECs Bq, Wq and (3 can be different in the valence and sea sectors. 
These equations tell us (to leading order in ChPT) what it means to set 
renormalized valence and sea quark masses equal: one chooses rriy^s such that 
Myy = Msg. We note that Wq^^s = if order-a improved Wilson fermions are 
used, and that Wov,s = as well as /3^,s = for GW fermions. 

Flavor-neutral propagators are given by {i and j are flavor indices) 
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with 

R=- (Mi - Ml) + (7,, + 7,, - 27,,) a' . (8) 

It is clear that R ^ even if M^^ = M^s- In fact, one may either choose the 
valence quark mass such that i? = 0, or choose it such that M^^ = Mgg and 
live with a non-vanishing R. Either way, the appearance of a double pole 
due to scaling violations in mixed actions is relevant for hadronic quantities 
sensitive to the double pole, especially if the effects are enhanced because of 
the IR-singular nature of the double pole. 



3 Example: 7 = pion scattering 



It is easy to find quantities particularly sensitive to the double pole. Exam- 
ples are the / = two-pion energy in a finite volume (which are related to 
the corresponding phase shifts) the anomalous large-time behavior of the 
Oo propagator [H], and the anomalous large-distance behavior of the nucleon- 
nucleon potential [7]. Here we will just briefly consider the finite- volume 
two-pion energy; for the Oq we refer to Ref. 

The energy of two pions at rest in an / = state, and in a finite spatial 
box of linear dimension L with periodic boundary conditions, is given by (to 
one loop) 13 12] 



2M„, 



where 



R 



Ml 



Sn^P ' 167r2/2 ' 

-0.53 + O(l/(ML)2) , 

49.59/(ML)2 + (l/(ML)3 



(9) 



(10) 



Setting M„ 
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Bo{ML) 
Ao{ML) 

■- Mss = M, we may consider two different regimes, one in 



which e ~ M^/A^ ~ uAqcd P!, and one in which e ~ M^/A^ ~ (aAQCDY 
jHj, where A ~ 1 GeV is the chiral symmetry breaking scale. In the first of 
these two regimes, the ratio of the one-loop to the tree-level terms in Eq. Q 
is of order x [ML)^ and x ML for the S"^, respectively, 6e terms. In the 
second regime this ratio is of order e x (ML)^. 
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One notes that the one-loop contributions are suppressed by positive pow- 
ers of e in both regimes, but for the positive powers of ML which accompany 
these powers of e. This is an example of the enhancement due to the "sick" IR 
behavior of the double pole — these terms would all disappear when i? = 0. 
But in a mixed- action theory at a 7^ 0, i? 7^ even at M^^ = Mgs- 

To get an idea about the size of these effects, one may substitute some 
typical values for the parameters of a LQCD simulation. Choosing for ex- 
ample aAqcD ~ 0.1, aM ~ 0.2 and L/a ~ 32, one finds that the S terms 
in Eq. (jH)) represent scaling violations of order 10%. This is small, but not 
negligible. Of course, in this estimate, we took all other constants (including 
the coefficient of in R) equal to one. Actual simulations are needed to 
make further progress with the investigation of the size of these effects. 

4 Conclusion 

Let us summarize what we learned about the use of mixed actions for LQCD 
computations. 

We begin by re-emphasizing that we assumed that the unphysical effects 
of mixed actions are encoded in scaling violations, and that EFT techniques 
can be used to study the issue. We recall that this approach has proven very 
successful in similar cases, such as the unphysical effects due to (partial) 
quenching. It is therefore important to test this assumption also in this case 
through actual simulations, in which in addition the numerical size of the 
effects can be estimated more reliably than through our parametric estimates 
above and in Ref. [2]. 

We have concentrated on the role of the double pole, because it is the 
most IR-sensitive probe of unphysical effects due to the use of mixed actions, 
just as it is in the case of (partial) quenching. Clearly, the effects are quantity 
dependent, and are expected to be most pronounced for those quantities for 
which there is a clear unphysical enhancement. The example we gave is that 
of enhanced finite-volume effects in pion scattering. We concluded in that 
case that the effects are likely to be rather small numerically in present-day 
simulations, but not so small as to be automatically negligible. 

Finally, we observe that in the case of the / = two-pion energy, this 
effect can also be monitored if one leaves out the "double-annihilation" dia- 
gram (as was done in Ref. .Q] because of poor statistics), because the effective 
theory can be adapted to reflect this situation For the terms shown in 
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Eq. this corresponds to dropping the 5"^ term. Our 10% estimate above 
came largely from the 5e term, and it is thus plausible that the size of the ef- 
fect does not become much smaller without the double-annihilation diagram. 
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